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I. INTRODUCTION 

There has been much interest recently in the phe- 
nomenon of dimerization in S' = 1/2 Heisenberg antifer- 
romagnets, where pairs of neighbouring spins couple to 
form 5 = singlet dimers. The dimerization may arise 
due to inhomogeneous bond interactions, as in the alter- 
nating Heisenberg chain (AHC) model, or the Shastry- 
Sutherland model in two dimensions jj]. Alternatively, 
it may emerge spontaneously, as the result of frustra- 
tion this seems to occur in the J1-J2 square lattice 
model at intermediate coupling values, for instance, al- 
though there is disagreement as to whether the pattern 
of dimerization is ordered ('valence bond solid') 3, 4] or 
disordered ('valence bond liquid' or 'resonating valence 
bond') HH. 

To understand the properties of dimerized phases, it 
is useful to construct an appropriate lattice formalism 
describing the dimers and their spin-triplet excitations. 
The physics of the system can then be connected with 
the properties of the elementary triplet excitations; and 
one can also use the formalism to construct a continuum 
'effective Lagrangian' field theory for the system at hand. 
Such a formalism was the 'bond-operator' representation 
constructed by Sachdev and Bhatt 7] (see also Chubukov 
0) some years ago, which is analogous to the spin- wave 
representation traditionally used to describe the magnet- 
ically ordered phases of these systems . 

Sachdev and Bhatt considered two spins Si and S2 
at either end of a single bond on the lattice, forming a 
dimer. They introduced a singlet and three triplet boson 
creation operators to form the corresponding states from 
the vacuum: 



s > 



\tx > 



\ty > 



st|o>^_(| Ti>-UT>) 
tl\o>=-l={\ TT>-Ui>) 
4|0>=-^(| TT>+IU>) 
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\tz> = 4lo>=^(i n>+UT>) (1) 

Then the spin operators can be represented 



S2a = ^[-sha - fls - ieapyt^pt^] (2) 



(where a, (3, 7 take values x, y or z), with the constraint 
that physical states must satisfy 



1. 



(3) 



They applied this formalism to develop a mean field the- 
ory of the frustrated square-lattice antiferromagnet. 

The problem with this approach is that the constraint 
Ol is awkward to implement analytically. Kotov et al. 
[lOj have applied an alternative "Brueckner approach" , 
in which the singlet operator is discarded, leaving only 
the constraint that two triplet excitations are not allowed 
on the same site (bond) . This is implemented by an infi- 
nite on-site repulsion term between triplets, which is ap- 
plied using an analytic Brueckner approach, valid when 
the density of triplets is small. The approach has been 
applied to the two-layer Heisenberg model [lol [Tlj , the 
quantum spin-ladder 12, 13], and the dimerized Heisen- 
berg chain with frustration [1J| , and some useful physical 
insights have been obtained. In particular, the occur- 
rence of two-particle bound states formed from the ele- 
mentary triplet excitations seems to be generic in these 
models. Nevertheless, the Brueckner implementation of 
the on-site repulsion term is also somewhat awkward to 
apply, and difficult to carry through in higher orders. 

Here we present an alternative approach in which the 
triplet exclusion constraint is implemented automatically 
by means of projection operators. We also use a "modi- 
fied" formalism, analogous to modified spin-wave theory 
[isinj^ . in which the two-body terms in the Hamiltonian 
are diagonalized through to the highest order calculated. 
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We denote the singlet ground state as |0 >, and introduce 
triplet creation operators that create the triplet states 
out of the vacuum |0 >, as follows 



FIG. 1: The alternating Heisenberg chain. 



The absence of any constraint makes the formalism eas- 
ier and more transparent to apply. The only drawback is 
the appearance of extra many-body interaction terms in 
the Hamiltonian, so that carrying the calculation to high 
orders would require the aid of a computer. 

To illustrate the formalism, we apply it to the case of 
the alternating Heisenberg chain (AHC). This model has 
itself attracted much attention recently, as new materials 
such as Cu(N03)2.2.5D20 [Hill have been constructed 
which appear to conform to this simple model, while at 
the same time more powerful neutron scattering facilities 
are coming on-line to explore their properties. For a re- 
view and further references, see Barnes et al. On 
the theoretical side, Uhrig and Schulz used a field 
theory approach to predict the appearance of both sin- 
glet {S = 0) and triplet {S — I) bound states below the 
two-triplet continuum. This was confirmed by later stud- 
ies [M, 21, 22]. Bouzerar and Sil 23] and Shevchenko et 
al. [ij] have treated the AHC using the Brueckner ap- 
proach; while Singh and Zheng |2^], Trebst et al. [23 
and Zheng et al. [Ig have carried out high-order dimer 
series expansions for the model, which give an accurate 
numerical picture of the dimerized phase. 

In Section II, we lay out the triplet-wave expansion for- 
malism for the case of the alternating chain. In Section 
HI, the expansion to leading orders on powers of the cou- 
pling A is discussed for the ground state energy and en- 
ergy gap. In Section IV, numerical results are presented 
for the ground-state energy, the one-particle spectrum, 
the two-triplon bound states, and the exclusive structure 
factors for these states. A summary and conclusions are 
presented in Section V. 
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Then the spin operators Si and S2 can be represented 
in terms of triplet operators by 



Sla 



r)[^a(-^ ^7^7) ^" (-^ tlty^ta i^afi^t\t. 



S20. = -[-ti{\ - t\t^) - {\ - t\t^)to, 



(7) 



where a, /9, 7 take the values x, y, z and repeated indices 
are summed over. This is similar to the representation 
of Sachdev and Bhatt 0, except that we have omitted 
singlet operators s\ s, but used projection operators (1 — 
tt^t^) instead. Assume the triplet operators obey bosonic 
commutation relations 



[tot : t^i^] — 5 



(8) 



then one can show that within the physical subspace (i.e. 
total number of triplet states is or 1), the representation 
obeys the correct spin operator algebra 

[Sia,Sifs] = iea/SfSij, [»S'2q,S'2^] = iea/37<S'27, (9) 



II. 



TRIPLET- WAVE EXPANSION 



The Hamiltonian for the alternating Heisenberg chain 
can be written 



= ^ Si • Si+i + A J2 Si • Si 



+1 



(4) 



i odd 



For A = 0, the system reduces to independent dimers as 
shown in Figure ^ Let us consider a single dimer with 
two spins Si,S2. The four states in the Hilbert space 
consist of a singlet and three triplet states with total 
spin S' = 0, 1 respectively, and eigenvalues 



Si • S2 — 



-3/4 (3 = 0) 
-1/4 {S = 1) 



(5) 



[Sla, S213] — 



(10) 



Sl = Sl = 3/4, Si • S2 = tlta - 3/4 (11) 

The projection operators ensure that we remain within 
the subspace. 

Returning to the alternating chain, we can now de- 
fine triplet operators tl^^ , tna for each dimer n along the 
chain. For a chain of N dimers, the Hamiltonian now can 
be expressed in terms of triplet operators as 
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n n 

n n 



This expression includes terms containing up to 6 triplet _ ^ J_ji/2 g-ifen^t ^]^3'j 

operators. For the purposes of the present calculations, " ^ ^ 

we shall drop terms with more than 4 triplet operators 
hencef or wards. 

Next, perform a Fourier transform 

^ (wc set the spacing between dimers d = 1), then the 

ifca = (^)^^^ Hamiltonian becomes 



fe fc ^'^''^ 123 

+47*la*2/3] + ^ X!'f'^l+2+3-4*L*L47*47 COS fci + ^1-2-3-447*27*30*40 COS ^4 

1234 

+(^1+2-3-4 [4a47*37*4a COS ^4 + 4a47*37*4a COSfci]} 

^Tm H ^1+2-3-4 4^4/3*3/3*47 - A^Afi^il^^ll COs(fcl - ^3) (14) 
1234 



Finally, as in a standard spin-wave analysis, we per- 
form a Bogoliubov transform 

*fca = CfeTfea -I- Sferl^.^ (15) 

where Cfe = cosh^^;, = sinh^^, 9-k = 0k, which pre- 
serves the boson commutation relations 

[TfeajTfe/^] = Skk'Saj3 (16) 

and is intended to diagonalize the Hamiltonian up to 

quadratic terms. After normal ordering, the transformed 
Hamiltonian up to fourth order terms reads 

H = Wo + H2 + H3 + Hi. (17) 

Here the constant term is 

Wo = 3Af[-^ + i?2 

-^(i?3 + i?4)(l - 8i?2 - 2i?i + i?3 - i?4)I18) 



expressed in terms of the momentum sums 

^2 = 

k 

k 

Ri = ^E^cosfc (19) 

k 

The quadratic terms are 

H2 = Y^i^krl^^ka + Qk{TkaT-ka + rLrlfc^)] (20) 

where 
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Ek = (c^ + Sfe)[l + 4A(i?3 + i?4)-2 cosfc(l-2i?i-8i?2-2i?4)]-AcfcSfc[cosA:(l-8i?2-2i?i + 2i?3)-2(i?3 + i?4)] (21) 

Qk = CfeSfc [1 - ^ (cos fc(l - 2i?i - 8i?2 - 2i?4) - 8(i?3 + i?4)) - ^ (4 + s^) [cos fc(l - 8i?2 - 2i?i + 2i?3) - 2(i?3 + Ri)] (22) 
The third and fourth-order terms are 



H3 = -^^eaf}^ ^['^1+2+3*3^^ (TlaT2^T3^ + 4^4p4a) + '5l+2-3*3^^ (T'L^p'^a'y + Tl'r'^2f3ria)] (23) 



and 



AN 



XI ['^1+2+3+4*4 (^L^la ^37^47 + ^Ia^2ar37r47) + (5i+2-3-4 (^f T|„r2^„T37T47 + ^4 A^T^^TaaTAj) 



1234 



^l+2+3-4«'4 (^LT2aT3V47 + Tl^T3^T2aTia) 



(24) 



where we have used the shorthand notation 1 • • • 4 for 
momenta ki ■ ■ ■ k^, and the vertex functions $3*' , $4*' are 
hsted in Appendix A. 

The condition that the off-diagonal quadratic terms 
vanish is 



This would leave some residual off-diagonal quadratic 
terms, arising from the normal-ordering of quartic op- 
erators. In a 'modified' approach M 

, we demand that 

these terms vanish entirely up to the order calculated, 
giving the modified condition 



(25) 



In a conventional spin-wave approach, this would be im- 
plemented in leading order only, giving the condition 



tanh 26k 



2ckSk Acosfc 
cl + si ~ 2[1 - A/2 cos fc] 



(26) 



tanh26k = 



A[cos fc(l - 8R2 - 2Ri + 2i?3) - 2(^3 + Ri)] 
2[1 - AcosA:(l - 2i?i - 8i?2 - 2i?4)/2 -f 4A(i?3 + i?4)] 

I 



(27) 



Self-consistent solutions for the N equations (|27(l . with 
the four parameters Ri ■ ■ ■ R4 given by equation p9|) , can 
easily be found by numerical means, starting from the 
conventional result H26() . 



III. EXPANSION IN POWERS OF A 

As a first check on the formalism, one may calculate the 
leading terms in an expansion of the energy eigenvalues 
in powers of A. From equation (|26|l . we easily see that to 
order A^ 

A A^ 9 

Sfe — — cos A; -I- — cos k 
4 8 



1 



A^ 
32 



cos^ k 



(28) 



and hence the lattice sums H19(l can be evaluated 

A^ A^ 
^16' 32' 

^ + 0(A^), R4 = 0{X') 



i?3 = 



(29) 



The leading-order behaviour of the vertex functions may 
easily be deduced from Appendix A. 

Substituting in equation H18|l , the ground state energy 
per site is 



Wo 
2N 



R2 



A 



(i?3+i?4)(l -8i?2 



'2Ri + i?3 — R4) 

3 3A^_ 3A3 
' 8 64 256 



A ^ (30) 



5 



1 /\3/\4 





FIG. 2: Perturbation diagrams contributing to the ground- 
state energy. 






FIG. 3: Perturbation diagrams contributing to the one- 
particle energy. 



in agreement with dimer series expansion results previ- 
ously obtained for this model One can easily show 
that perturbation diagrams such as those in Figure |2] do 
not contribute until O(A^) or higher. 

The energy gap at leading order can be found from 
equation J^IJ: 



Eh 



— cos k + — [i — cos^ k] , A 
2 8 







(31) 



The perturbation diagrams Figures Ot) and l^ls) also 
contribute at order A^. Note that diagram (SJa) does not 
appear in the formalism of Shevchenko et al. |14| : the 
extra terms in our formalism are needed to implement 
the hardcore constraint that two triplons cannot occupy 
the same site. At leading order, the contributions of these 
diagrams are 



AE 



3a) 



A2 



(1 + cosfc), A 







(32) 



which again agrees with series expansion results |24j . 

If we compare equation 1)34(1 at small momentum with 
the continuum dispersion relation for a free boson, 



Efc ~ \/ rrfid^ + k^c? 



(35) 



we readily discover the leading behaviour of the effective 
triplon parameters, i.e. the triplon mass 



TO - -[1 + A + 0(A2)] 
A 



and the 'speed of light' 



c-^-^ + 0(A3) 



(36) 



(37) 



in lattice units. Note that the mass diverges and the 
speed of light vanishes as A ^ 0. 



IV. NUMERICAL RESULTS 



Writing the Hamiltonian as 



where 



and 



H = Ho + V 



(38) 



(39) 



(40) 



we can treat Hq as the unperturbed Hamiltonian and V 
as a perturbation to obtain the leading-order corrections 
to the predictions for physical quantities outlined in the 
previous section. Numerical results for the model have 
been obtained using the finite-lattice method. The mo- 
mentum sums are carried out for a fixed number of dimers 
A'^, using corresponding discrete values for the momen- 
tum fc, e.g. 



kn — 



2Tm 



n = l.---N 



(41) 



Results were obtained for N up to 40, and a fit in powers 
of 1/iV was made to extrapolate to the bulk limit N 



AE] 



3b) 



A^ 
' 4 ' 



A^O 



(33) 



(see the next section for further details). This gives a 
total single-particle energy 



1 



A 



■ cos k 



A2 



cosfc[2-t-cosA:], A -> (34) 



A. Ground-state energy 

The leading corrections to the ground-state energy cor- 
respond to the diagrams in Figures and [2b) • Their 
contributions are 
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Ae 



2a) 



123 



$£;(123)$ni23) 
{El +E2+ Es) 



-51 



+2+3- 



(42) 



Ae, 



26) 



-3A2 

47V3 



1234 



1+2+3+4 



$^^(1234) 



(-El +E2 + E3 + Ei) 



3$i^^(1234) + $^'^(1324) + $^'^(1423) 



.(1) 



(1)^ 



(43) 



In leading order one can show that these terms are 0(A'*), 
whereas diagrams such as FigurelSt) are O(A^) or higher. 
The resuhing bulk estimates of the ground-state energy, 
including these corrections, are listed in Tabled Figure 
2] shows the behaviour of the ground-state energy as a 
function of A resulting from this modified triplon theory, 
as compared with the high-order dimer series calculations 
of Zheng et al . l24l and exact diagonalization data of 
Barnes et al. [l^l- It can be seen that out to A ~ 0.4 
there is quantitative agreement between our calculation 



and the series estimates, but some discrepancy emerges 
at larger A. 



B. One-particle spectrum 

The leading corrections to the one-particle spectrum 
correspond to the diagrams in Figures|2tt) andl^J)). Their 
contributions are 



12 



(Ek -El- E2) 



(44) 



Ai?f) 



$f^(123fc) 



El — E2 — E3 



3$i'^'(123A:) + $^"^(321fc) + $^"'(312fc) 



(4)/ 



(45) 



In leading order, these terms are 0{X^), as stated in the 
previous section, while diagrams like|3t), d) are 0{\'^) or 
higher. 



The resulting bulk estimates of the energy gap at /c = 
are listed in Tabled and displayed in Figure |S1 It can 
be seen that the inclusion of the diagrams |2K) and OJd) 
improves the agreement with series dramatically. This 
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TABLE I: Values for the energy per dimer eo and the energy 
gap at A; = as functions of A. The left-hand box giver 
series estimates (2^ while the right-hand box gives our present 
trip let- wave results. 



1.15 



A 


Series 


expansion 


Triplet 


expansion 








eo 


Energy gap 




eo 


Energy gap 








-0 


75000 


1.00000 


-0 


75000 


1.00000 





1 


-0 


75096 


0.94628 


-0 


75096 


0.94647 





2 


-0 


75394 


0.88521 


-0 


75392 


0.88625 





3 


-0 


75914 


0.81684 


-0 


75896 


0.81885 





4 


-0 


76672 


0.74106 


-0 


76611 


0.74252 





5 


-0 


77694 


0.65748 


-0 


77535 


0.65483 





6 


-0 


79010 


0.56530 


-0 


78659 


0.55341 





7 


-0 


80662 


0.46300 


-0 


79970 


0.43649 





8 


-0 


82712 


0.34753 


-0 


81455 


0.30312 





9 


-0 


85268 


0.21130 


-0 


83096 


0.15309 


1 





-0 


88630 


0.00828 


-0 


84878 


-0.01323 




FIG. 5: [Color online] Energy gap at A; = as a function of A. 
The dot-dashed line shows the series estimates [2^ . while the 
other lines show the leading order and improved triplet-wave 
results. The filled squares are results from Shevchenko et al. 
[T^. and circles are results from Barnes et al. |Tg|. 




(a) (b) 

FIG. 6: One-particle spectrum at selected couplings. 
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FIG. 7: a) The inverse triplon mass parameter 1/m, and b) 
the 'speed of light' c^, as functions of A. 



note that the exact value of the speed of light c at A = 1 
is 7r/2 = 1.57 [jj; which is about twice the value of even 
the series estimate (~ 0.78). This is presumably due to 
the singular behaviour of the model in this limit, includ- 
ing logarithmic corrections, which even high-order series 
expansions cannot accurately reproduce. 



agreement may be fortuitous, given that the agreement 
for the ground-state energy is not so good, but it is grat- 
ifying to see nevertheless. It can be seen that our present 
approach improves upon that of Shevchenko et al. [l^ 
at intermediate A. 

The dispersion of the one-particle energy as a function 
of momentum is illustrated at selected couplings in Fig- 
ure while Figures |7^) and [3d) show the corresponding 
behaviour of the inverse triplon mass parameter 1 /m and 
the speed of light squared, c^. At the smaller coupling, 
the dispersion agrees quantitatively with series estimates, 
but at A = 0.8 we can see that the minimum of the en- 
ergy is too broad: the curvature at fc = should diverge 
as A — > 1. This is reflected in the fact that our results 
for 1/m and are much too low at large couplings. We 



C. Two-triplon bound states 

It has been found in previous studies [13. that the 
quartic terms in the Hamiltonian lead to attraction be- 
tween two elementary triplons, giving rise to S* = and 
5 = 1 bound states. We look for solutions of the two- 
body Schrodinger equation 

H\^ >= E\^j > . (46) 

The two-body wave functions \^p{K) > can be written 
as follows: 

Singlet sector (S = 0): 



8 



(47) 

where K is the centre-of-mass momentum and q the rel- 
ative momentum of the two particles; 

Triplet sector (S = \): 



wave functions: 



l^A/^^.^(JC,9,p)V'5'^(X,p). (49) 



K/2+q,f3'^K/2^q,y\^ > 



(48) 

where K is the centre-of-mass momentum and q the rel- 
ative momentum (we will not write out the quintuplet 
states explicitly). 

From equation (|46l) one can readily derive the inte- 
gral Bethc-Salpeter equation satisfied by the bound-state 



In leading order, the scattering amplitudes 
M^''^{K,q,p) are simply given by the 4-particle 
vertex from the perturbation operator V, Figure |SJi)- 
Hence we find for the different sectors: 



Singlet sector (S — 0): 



{K,q,p)^- [3$f ' {K/2 + p, K/2 -p,K/2 + q, K/2 



q) + 'i>'i^^iK/2+p,K/2-p,K/2 + q,K/2~q)] (50) 



where the wave function is symmetric, 

t^^iK,-q)^i;''{K,q) (51) 
and the symmetric and antisymmetric pieces of the ver- 

(3) 

tex function $4 are defined: 

= ^(1234) ± $f\l243)] (52) 
Triplet sector (S — I): 



and 



A 



AP\K,q,p) = iK/2+p,K/2-p,K/2+q,K/2-q) 

(53) 

with the wave function antisymmetric 

^^iK,^q)^-i;'^iK,q). (54) 
Quintuplet sector (S = 2): 

M^{K,q,p) = ^$f^+(i\:/2+q,ii:/2-q,ii:/2+p,i\:/2-p) 

(55) 

where the wave function is once again symmetric 

^'^{K,-q)^i;'^{K,q). (56) 

At leading order in A, we find 

{K, q,p) ^ A[cos(iir/2)(cosp + cos 5) — cos p cos g] 

(57) 



M'^{K, q,p) ^ — — sinpsing 



(58) 



Following Shevchenko et al. , one can then find sim- 
ple solutions (unnormalized) to the Schrodingcr equation 
(El: 



5 cos(iir/2) - cosq 

^ ^ 1 + cos2 {K/2) -2 cos( A72) cos q ' 



i^^{K,q) 



sm q 



,(59) 



1 4cos2(i^/2) - 4cos(is:/2) cosg 
corresponding to bound-state energies 

E^{K) - 2- ^(l + cos2(i^/2)) 

E'^iK) - 2- ^(l + 4cos2(ii'/2)) (60) 

compared to the lower edge of the 2-particle continuum 
E2{K) ^2-\cob{K/2). (61) 



These agree with the dimer series expansion [I3, |23 at 
leading order. 

Note, however, that the singlet solution is valid for all 
K ^ touching the continuum at if = 0; while the triplet 
solution is only physically valid over a finite range of 
momenta K = {27r/3, 47r/3}, and not at K = 0. The 
end-points of this range are just where the triplet bound 
state enters the continuum, and the denominator of H59|) 
for ip'^ vanishes at q = 0. Note also that both dispersion 
curves meet the lower edge of the continuum at a tangent. 



9 





/ \ 



(a) 



(b) 
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FIG. 8: Perturbation diagrams contributing to the 2-particle 
scattering amplitude. 



At the next order O(A^), further diagrams contribute, 
as given in Figures Two of these, Figures|Hh) and 
IHJ), we are not in a position to calculate, because they 
involve 5 or 6-particle vertices. Figure |H|3) is already ac- 
counted for by diagonalizing the effective Hamiltonian 
in the 2-particle subspace. Diagrams |Hli) and sim- 
ply correspond to renormalizations of the single-particle 
energies in the diagonal terms of the effective Hamilto- 
nian in the 2-body sector. Finally, we can calculate the 
contribution of Figures IS]:), |Ht) and|S^) to the effective 
Hamiltonian using perturbation theory. In general, the 
change in the energy eigenvalue is 

^E=^Y1 ^^^(^' q,PWK,PmK, q) (62) 
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where the vertex function AM{K,q,p) for each different 



FIG. 9; Dispersion relation for the triplet {S — 1) 2-particle 
bound state at selected couplings a) A = 0.05 and b) A = 0.25. 
The lower edge of the continuum is labelled 'cont'. 



diagram and spin state is listed in Appendix B. 

The corrections due to these diagrams can now be cal- 
culated. On a finite lattice, equation (|49l) becomes a 
matrix eigenvalue equation, which can readily be solved 
numerically. We have calculated results for lattices of 
up to = 40. The resulting bound-state spectrum is 
displayed in Figures El and ^| The first thing to note 
is that the modified but uncorrected triplet expansion 
agrees with series expansion estimates quite well, for both 
the lowest-lying singlet and triplet bound states. For the 
triplet state, the result is substantially better than that 
of Shevchenko et al. 0| at A = 0.25. Inclusion of the 
perturbation corrections actually makes the agreement 
worse, and gives much too large binding energies, espe- 
cially at A = 0.25. This can be attributed to the neglect 



10 



2.01 
2 

1.99 
1.98 
1.97 
1.96 
1.95 
1.94 



1.95 
1.9 
1.85 



Series expansion 
Triplet expansion 
Triplet exp. w/ corrections 



cont. 



S=0 



0.5 



1 1.5 
K 

(a) 



2.5 



continuum 



Series expansion 
Triplet expansion 
Triplet exp. w/ corrections 
Shevchenko et al 




(b) 



FIG. 10: [Color online] Dispersion relation for singlet {S = 0) 
2-particle bound states at selected couplings a) A = 0.05 and 
b) A = 0.25. The lower edge of the continuum is labelled 
'cont'. 



of diagrams which are of the same order as the 

diagrams we have calculated. Unless the extra diagrams 
are included, we cannot do better than the uncorrected 
estimates. 

We have also looked for signs of the second singlet and 
second triplet bound states which were found to appear 
at order by Trebst et al. In the corrected results, 
a second singlet bound state does appear, in fact, but 
with much too large a binding energy once again. The 
detailed dynamics of the bound states are sensitive to 
higher-order terms. 

D. Structure Factors 



The "reduced exclusive structure factor" or spectral 
weight for a specific intermediate state A with momen- 
tum K can be written 



(63) 



where 



(64) 

and the sum i* runs over sites of the unit cell on the 
lattice, and N is the number of unit cells (dimers) . Using 
equations Q, H13|) . and IjlSfl . the spin operators 5*1 and 
^2 on sites 1 and 2 can be expressed in terms of triplet 
operators (taking n = in equation Ijl^fl ): 



k 12 

[Tt\mirlrlTl + T^rln^ + tIt^^t^^ + t^.t^.t^^) + (123) (rt^rt^Ts^ + TlT2^n^)\{(y^) 



123 

r 



where the upper and lower signs correspond to 81,82 r2^''(12) = — — (ciS2 — S1C2) (67) 

respectively, and 87V 



T\'\k) = -^(cfc + sfe)(l - i?i - AR2) (66) 



2ViV 
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FIG. 11: Perturbation diagrams contributing to exclusive 
structure factors. 



. , Ka . , A 
zsml ) 1 H — cos A 



'64 



(5cos2X-7)]. 



(71) 



Here a represents the spacing between spins in the dimer, 
i.e. a — 1/2 for the uniform lattice in our present units. 

Higher-order diagrams such as Figs. Illb). c) do not 
contribute until O(A^). Their contributions are listed in 
Appendix C. Hence we find 



(68) 



^tp^^^ - zcos(^)y sinJs: 



(72) 



r3^^^(i23) 



(ci + Si) , 

4iV3/2 



(,C2S3 + S2C3J 



(69) 



and 



(70) 



In leading order fFigurelTTk'). the one-particle matrix 
element is 

Ka 

ni{K) = ismi — )[il-R,-AR2){cK + SK)] 



^a(llc) 

Ip 



I sm ) — . 



(73) 



We must also account for the renormalization of the 1- 
particle wave function due to Figures|3Ji) andl^b), giving 
a multiplicative renormalization factor 



Z 



K 



2N 



1+2-K 



12 

+$i^)(321X) + $i'*)(312X)] 

1 cos K 

8 16 



Ek-Ei- E2 



16N' 



123 



1+2+3-K 



$i^'(123if) 



{Ek -E1-E2- Esf 



giving a total amplitude 



,Ka 



^%{K) - i[sin( — )[l-f jcosi^-l-— (-ll-8cosis: + 5cos2ii')]+cos( — )— siniiT] 



whereas Zheng et al. |2g obtain 



,Ka 



n^piK) - i[sin( )[l-f -cosi^-l-— (-ll-4cosi4'-|-5cos2i4')]-|-cos( )— siniiT] 



(74) 



(75) 



(76) 



We have been unable to resolve the source of the discrep- 
ancy at order A^. 

The calculated results for the one-particle spectral 
weight are displayed in Figure 12. It can be seen that 
the results match the series estimates quite well, even at 
larger A. The corrected estimates are a little low at small 



k, and a little high at large k: this reflects the discrep- 
ancy at order A^ referred to above. 

For the 2-particlc bound states, the leading perturba- 
tion diagrams contributing to the exclusive structure fac- 
tors are illustrated in Figs. IlltD.e). Their contributions 
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for the triplet states are 



n'^^^^'^^K) = -8^ViVcos(^) X 



J2 ^'^(K, qM'^ iK/2 + q, K/2 - q) 



-I — — cos( — ) sm{K 2) 
2VN ^ 2 ' ' ' 

-."^ sin qi;^{K, q) 



(77) 



at the threshold points {cos(iir/2) = ±1/2} where the 
bound state merges with the continuum. 

The numerically calculated results for the spectral 
weights are displayed in Figure 12. For the 1-particle 
weight, it can be seen that the corrected triplet-wave ex- 
pansion matches the series estimates very well at the two 
lower couplings, and only begins to deviate significantly 
at A = 0.5. The triplet wave expansion also works sur- 
prisingly well for the 2-particle weight, which is only of 
order 1% of the 1-particle weight. 



and 



= 2iAsin( — )Ti^'^(i^) 



K/2+q + ^ 



A ,Ka^ , , , 
i^smf— )cos ii:/2) 



^sinqiFiK, q) 



(78) 



Inserting the wave function EHl we obtain the leading 
order behaviour of the triplet bound state contribution 
to the structure factor (for the uniform case a — 1/2): 

\2 

S^'^iK) - y sinS(iC/4)(l-4cos2(i^/2)), (79) 

which agrees with the leading order series calculation 
[26| . Thus the bound-state spectral weight vanishes 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have developed a modified triplet- 
wave expansion method for dimerized spin systems, anal- 
ogous to the modified spin-wave formalism 0, 0| for 
magnetically ordered systems. It differs from the ear- 
lier app roaches of Sachdev and Bhatt 7] and Kotov et 
al. [lO| in that projection operators are used to confine 
the system to the physical subspace in the bosonic for- 
mulation, eliminating the need for a separate constraint. 
The two-body boson operators are also fully diagonalized 
through the highest order calculated. 

The formalism has been applied to the case of the alter- 
nating Heisenberg spin chain. Using perturbation theory 
to second order, we have calculated the ground-state en- 
ergy per dimer, the dispersion relations for one-particle 
states and two-particle bound states, and the spectral 
weights for these states. It has been shown that the re- 
sults reproduce the leading order terms in a dimer series 
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expansion in powers of A p^. I25I , apart from an unex- 
plained discrepancy at order in the 1-particle spectral 
weight. The results are quantitatively accurate at small 
A, but begin to show significant discrepancies from high- 
order series expansions at larger A, as one would expect. 
The discrepancies become more serious for the more sen- 
sitive dynamical quantities such as two-particle binding 
energies. The inclusion of a partial set of higher-order 
corrections for the two-particle binding energies made 
things worse rather than better, as one perhaps should 
have expected: all terms of a similar order in A must be 
included simultaneously if a good result is to be obtained. 
Nevertheless, the qualitative behaviour is correctly repro- 
duced by the formalism. In particular, the formation of 
two-triplon bound states near K = 7r/2 in both the sin- 
glet {S = 0) and triplet (S = 1) channels, which was 
discovered previously [l^ l20l | is reproduced. 

The behaviour of the triplet bound state near the 
threshold where it merges with the continuum is interest- 
ing. We have seen that the bound-state dispersion curve 
merges at a tangent to the continuum, and that the spec- 
tral weight vanishes at the threshold. The bound-state 
solution does not extend into the continuum, but termi- 
nates at the threshold. This provides a neat example of 
the phenomenon of "quasiparticle breakdown" discussed 
recently in the literature '28, '29II33: i.e. the termination 
of a single-particle state where it enters the continuum 



for one-dimensional systems. 

Our results appear to be more accurate and reliable at 
intermediate couplings A than those of Shevchenko et al. 
jT^. However, they cannot match the quantitative accu- 
racy of the high-order dimer series expansions p^l25ll2^ 
or exact diagonalization on large lattices . The calcu- 
lations could be pushed to higher orders with the aid of a 
computer, but it is doubtful whether this is worthwhile. 
The main value of a 'lattice bosonization' approach such 
as this is to provide a better analytic understanding of 
the behaviour of the model, and a half-way house towards 
a continuum 'effective field theory' for the model. For in- 
stance, we have shown how the triplon mass parameter 
and the 'speed of light' can be calculated, which would be 
fundamental parameters of the effective field theory. It 
would be interesting to apply the approach to dimerized 
models in two dimensions. 
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APPENDIX A 



The vertex functions $3*'', $4*'' are: 



*3^^(123) = ^[sinfci(ci -f Si)(c2S3 - S2C3) + sinfc2(c2 + S2)(c3Si - S3C1) 

-|-Sinfc3(c3 -f S3)(ciS2 - S1C2)] (80) 

(2) 1 

$^'(123) = -[sinfci(ci + si)(c2C3 - S2S3) + sinfc2(c2 + S2)(siS3 - C1C3) -I- sinfc3(c3 + S3)(siC2 - C1S2)] (81) 

$4^''(1234) = i[(c0sfci + COsfc2)(ciC2 -f C1S2 + S1C2 + SiS2)(c3S4 + S3C4) 
+ (C0S fcs + C0SA;4)(C3C4 + C3S4 + S3C4 + S3S4)(ciS2 + S1C2) 
-I- COs(fci + fc4))(ciS4 — SiC4)(s2C3 — C2S3) 

-f COs(fci + fc3)(ciS3 - S1C3)(S2C4 - C2S4)] (82) 
(2) 1 

$4 (1234) = -[(cosfci + C0S/S2)(C1C2 -I- C1S2 + S1C2 -I- SiS2)(s3C4 -f C3S4) 

-|-(C0S fc4 -I- C0sA;3)(ciS2 + SiC2)(c3C4 -|- C3S4 + S3C4 + S3S4) 
-|-COs(fci - fc3)(siC2S3C4 -I- C1S2C3S4 - C1C2C3C4 - S1S2S3S4) 

-I- COs(fci - fc4)(ciS2S3C4 -|- S1C2C3S4 - C1C2C3C4 - S1S2S3S4)] (83) 
$4 •^(1234) = COsA:i((ciC2 + SiC2)(s3C4 + C3C4) -t- (ciS2 + S1S2)(S3S4 -I- C3S4)) + COsfc2((ciC2 -|- CiS2)(c3S4 + C3C4) 
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+ (siC2 + S1S2)(S3S4 + S3C4)) + COSfc3((ciS2 + SlS2)(c3S4 + S3S4) + (siC2 + CiC2)(s3C4 + C3C4)) 
+ COSfc4((siC2 + S1S2)(S3C4 + S3S4) + (ciS2 + CiC2)(c3S4 + C3C4)) 
+ COs(fci — fc4)(ciC4 — S1S4)(C2C3 — S2S3) 

+ COs(fci + fc2)(ciS2 - S1C2)(C3S4 - S3C4) (84) 

$4'^'' (1234) = (cosfci + COsfc2)(ciC2 + S1S2 + C1S2 + S1C2)(C3C4 + S3S4) + (C0sfc3 + COSfc4)(ciS2 + SlC2)(c3S4 
+ S3C4 + C3C4 + S3S4) + COs(fc2 + fc3)(ciS2C3C4 + S1C2S3S4 - C1C2S3C4 - S1S2C3S4) 

+ COs(fci + fc3)(siC2C3C4 + C1S2S3S4 - C1C2S3C4 - S1S2C3S4) (85) 

We have 'symmetrized' these expressions with respect to their indices, using momentum conservation. 

APPENDIX B 

The two-body terms AM{K,p, q) defined in equation H52I) for the diagrams Figs. |Ht), f) and g) are as follows (the 
energy denominators are 'symmetrized' between initial and final states): 

Scalar state 

AA4'\K,q,p) = (86) 
' <i>f (p - q, K/2 + q, K/2 + p)<f f {K/2 -p,p-q, K/2 - q) 



AMf\K,q,p) = \' 



{Ep^q + l/2(£';^;/2+9 + Ek/2-p - Ek/2+p - EK/2-q)) 

+ {ditto,p —p) + {ditto, q —q) + {ditto, p ^ —p, q ^ —q)} (87) 

AM(^^' {K,q,p) = - E { ^^i^ {p~q~k,k,K/2 + q,K/2+p) (3a>i^) {k,p~q-k,K/2-p,K/2-q) 

k 

+^f{K/2 ~p,k,p~q~k,K/2-q) + $f ^(if/2 -p,p-q-k,k,K/2- q)) 
+2<^^1\k/2 + q,p-q-k,k,K/2+p)x {3<i>[^\K/2 - p, k,p - q - k, K/2 - q) 

{k,p-q~k,K/2-p,K/2-q) + $f {K/2 ~p,p-q-k,k,K/2-q))] 
/ {Ek + Ep^q^k + 1/2(£'a:/2+9 + Ek/2-p — Ek/2+p — EK/2-q)) 

+{ditto,p <-> —p) + {ditto, q ^ —q) + {ditto, p ^ ~p, q ^ ~q)} 



Triplet state 



AMr\K,q,p) . ,2 f [K/2^p, K/2 - p, K)^^^^ {K/2 + K/2 - q, K) ^^^^ 

{Ek - l/2(£'K/2+g + EK/2-q + Ek/2+p + Ek/2-p)) 

rm),^ ^ _ A2 / <i>f ^ {p - q, K/2 + q, K/2 + p)<i>f {K/2 ^p,p^q, K/2 ^ q) 



AM^''>{K,q,p) = - 



2 1 {Ep^q + l/2{EK/2+q + Ek/2-p — Ek/2+p — EK/2-q)) 

{ditto, p ^ ~p) — {ditto, q ^ —q) + {ditto, p ^ —p, q ^ —q)} (90) 



AM^''^ {K,q,p) = - E { {p-q-k,k, K/2 + q, K/2 + p)(3<i>f {k,p-q^k, K/2 - p, K/2 - q) 

k 

{k,K/2~p,p-q-k,K/2-q) + <^'^^\p-q-k, K/2 -p,k,K/2-q)) 
+2^f' {p-q-k, K/2 + q,k, K/2 + p) x {^[^^ {k,p-q-k, K/2 - p, K/2 - q) 
{p-q-k, K/2 - p, k, K/2 - q))] 

/ {Ek + Ep^q^k + l/2(i?K/2+g + Ek/2-p - Ek/2+p - EK/2-q)) 

— {ditto, p —p) — {ditto, q ^ —q) + {ditto, p ^ —p, q ^ —q)} (91) 
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APPENDIX C 

Contributions to the 1-particlc matrix elements ^l'^{K) from the diagrams shown in Figure [TTb).c) are: 

<"l(i^) = A^X cos{^) Y: j^J'^^^l^l [^'i^ ^K) - $f (12 i^)] (92) 
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